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Abstract 

Electro-magnetic probes such as dileptons and photons are strong probes to investigate the thermodynamical state 
of the early stages of collisions since they leave the system unscathed. The PHENIX experiment has measured both 
photons and dileptons in p+p, d+Au and Au+Au collisions. An excess of dilepton yield over the expected hadronic 
contribution is seen in 0.2-0.8 GeV/c^ in Au-hAu collisions, which is prominent in lower pr and most central. Direct 
photons are measured through their internal conversion to electron pairs. We saw a large enhancement in Au+Au 
collisions over p+p yield scaled by the number of binary collisions. It turned out from the latest results on d+Au 
collisions that this enhancement is not explainable by a nuclear effect. 
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1. Introduction 

Many intriguing phenomena have been observed at RHIC where a hot and dense matter is expected to be created. 
The large suppression of the yield of the single hadrons at high transverse momentum {pj) suggested that the matter is 
so dense to stem fast partons with large emerged from the very early stage in the medium [1]. The large elliptic flow 
of particles and its scaling in terms of kinetic energy of the particles suggests that the system is locally in equilibrium 
at as early as 0.3 fm/c. Although the hadronic probes already exhibited many interesting phenomena, since they are 
suffered from strong interactions in later stages to some extent, observation of more direct and penetrating probes have 
been desired. 

Electro-magnetic probes such as photons or dileptons are ideal in this sense, since they interact with medium or 
other particles only electro-magnetically, once produced [2]. Therefore, these probes are of interest already from the 
beginning of the history of the heavy ion collisions. At the leading order, the production processes of photons are 
annihilation {qq — > yg) and Compton scattering {qg yq). Their yields are proportional to aas, which are ~40 times 
lower than the ones from strong interaction. The dilepton production is mainly from annihilation of quark and anti- 
quark {qq 7*), and the yield is proportional to a^. Except for some difference in production processes, photons and 
dileptons can be treated as same observables. A theoretical calculation tells that photons and dileptons share similar 
emission sources [3]. Fig 1 shows a mapping of photons and dileptons in one coordinate. 

PHENIX has measured dileptons and photons since the RHIC has started its running in 2000. In this paper, we 
show the most recent results on dileptons and photons in p+p, d+Au and Au+Au collisions at ^snn=200 GeV. 

2. Dilepton analysis 

The PHENIX detector has a capability of measuring momentum of charged particles with high accuracy and of 
identifying electrons with high efficiency and strong rejection power of hadrons [4]. In measuring dileptons, we have 
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Figure 1 : Photons and dileptons in one coordinate with three degrees of freedom. 



a huge combinatorial background e^e~ arising from random combination of electrons from photons converted at beam 
pipe, Dalitz decays of or r]. We constructed combinatorial background by combining electrons or positrons from 
different events, and subtracted them from foreground mass distributions. 

The subtracted distributions still include unphysical correlated pairs from back- to-back jets. We estimated the 
contribution using PYTHIA event generator followed by a detector simulation, and subtracted it off from the invariant 
mass distribution. The signal to background ratio is ~ 10"^ in total. The resulting distribution is corrected for 
efficiencies of single electrons, and compared to various hadronic contributions. 

3. Low mass and low pr dileptons in Au+Au collisions 

Fig 2(a) shows the invariant mass distribution of electron pairs measured in Au+Au collisions at y5Aw=200GeV [5]. 
The data is compared with electron pairs from various hardonic sources calculated using a Monte Carlo and filtered 
by the PHENIX acceptance. A large excess is seen in 0.2-0.8 GeV/c^ in Minimum bias Au+Au collisions. In order to 
study this excess in detail, we divided the data set into various pr, mass and centrality bins. Fig. 2(b) and (c) show the 
mass spectra in various pr and centrality bins, respectively. The excess dies out at higher pr and lower centralities, 
which suggests that the excess is contributed from a thermal source. We made summary plots for the detail studies as 
shown in Fig. 3(a) and (b). As we already saw in the invariant mass distributions, the yield of excess region increases 
as collisions become more central. The slope parameters are low for low mr - mo and high for high - mo, implying 
that the local slopes are roughly proportional to the energy of photons emitted. This is another evidence of that they 
are from thermal sources. 



4. Low mass and high pr dileptons in Au+Au and p+p 

Turning the region of interest to higher pr, where pr » M, the yield of dileptons are considered to be dominated 
by internal conversion of real photons as depicted in Fig 4(a). Taking this advantage, we performed direct photon 
measurement through dilepton measurement [6]. The invariant mass distribution of Dalitz decay of r] and direct 
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Figure 2: (a, left) Invariant mass distributions in Minimum bias Au+Au events with hadronic cocktail calculation, and the ones in (b, middle) 
various pj regions and in (c, right) centrality classes. 
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Figure 3: (a, left) Yield and (b, right) slope parameters of dileptons in low mass and low pj region. For slope parameters, solid points are calculated 
from truncated means of the mass distributions, while open points are from exponential fits to the distributions. 




Figure 4: (a, left ) Diagram of real photons and internal conversion of photons, (b, right) Fit to data with direct photon and hadronic cocktail 
distributions. 
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photons can be calculated using Kroll-Wada formula as shown below. 
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Then, the real invariant mass distribution is fitted with a function of: 

F = (1 - r)f, + rfd 

where fc is the cocktail calculation, fd is the mass distribution for direct photons, and r is the free parameter in the fit. 
Using Kroll-Wada formula, r is converted to direct photon to inclusive photon ratio as follows: 

JmA^ee > 0.15) y*diAniee ~ 0) y^-^ 

r = oc — — 



7* (mee > 0.15) rlMee ~ 0) 



inc 



Once is obtained, the invariant yield of direct photons are calculated as jinc x r^. The procedure demonstrated in 
Fig 4(b) is for 1.0< pj <1.5GeV/c. The dotted lines show the contributions from various hadrons, the solid blue 
is the sum of these contributions, and the solid red line show the distribution by direct photons converted internally. 
The r value is determined by the fit to the data. The error of the fit corresponds to the statistical error. We applied 
the procedure for various pt^ and centralities in p+p and Au+Au collisions, and obtained the pr spectra as shown in 
Fig 5(a). The result is obtained from 810 million Minimum bias Au+Au events collected in RHIC Year-4 run, and 




Figure 5: (a, left) Direct photon yields in p+p and Au+Au collisions, (b, right) Comparison of the yield in 0-20 % Au+Au collisions with various 
theoretical models. 

1.5 billion minimum bias events and 270 million high pj single electron trigger events collected in RHIC Year-5 p+p 
run. The distributions are for 0-20 %, 20-40 % centrality and Minimum Bias collisions in case of Au+Au collisions. 
The distributions were then fitted with the p+p fit plus exponential function. We obtain slopes and dN/dy for three 
centralities as shown in Table 4. There are many models in the market that have a wide range of initial temperatures 
and thermalization times. Fig 5(b) shows the comparison of the data with the models. Since the initial temperature 
and thermalization time are highly correlated, the comparison of data and models does not help constraining two 
parameters to definite numbers. 
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Table 1: dN/dy and slopes of direct photon pj distributions obtained from the fit to the data with the p+p fit plus exponential function. 



Centrality 


dN/dy {pT >lGeV/c) 


Slope (MeV) 




0-20 % 


1.50+0.23+0.35 


221+19+19 




20-40 % 


0.65+0.08+0.15 


217+18+16 


5.0/4 


MinBias 


0.49+0.05+0.11 


233+14+19 


3.2/4 



5. Low mass and high pt dileptons in d+Au 

It has been a question whether or not the excess is purely from the source that only exists in Au+Au collisions. A 
possible effect that increases the the yield in nucleus collisions is kr broadening, or so-called Cronin effect. In order 
to quantify the effect, we analyzed d+Au data from RHIC Year-8 run with the same procedure. Fig. 6(a) shows the 
invariant mass distribution in d+Au collisions for the lowest pr where the direct photon yield is significant. After 
repeating over pj bins, we obtained the Vy as a function of pr. Fig 6(b) shows the Vy in various collision systems. The 
lines show the ratios calculated using the direct photon contributions estimated by NLO pQCD prediction scaled by 
the number of binary collisions (Ncoii)- Three lines show the predictions with three mass scales (p = pt/2,pt, 2pt). 
The p+p ratio is very close to the NLO pQCD prediction, while Au+Au data largely deviates in low pr from the 
prediction. The d+Au data shows a moderate excess, but the excess is not as much as the one in Au+Au collisions. 
From the ratios obtained, we calculated the direct photon spectra, again following the procedure explained previously. 




Fig 7(a) shows the direct photon spectra for d+Au Minimum bias events, together with the real photon result from 
RHIC Year-3 run. A comparison is made with direct photon spectra from p+p collisions scaled by Ncoii- Fig. 7(b) 
shows a comparison of data with NLO pQCD calculation scaled by Ncou- In both cases, the yield in d+Au collisions 
are higher than the one expected from p+p yield, suggesting that the nuclear effect is seen in d+Au collisions. 

The ultimate interest in the series of the direct photon measurements is whether or not the excess seen in Au+Au 
collisions can be explained by cold nuclear effect (Cronin effect). In order to evaluate this, we compared the Au+Au 
yield with the d+Au yield scaled by the difference of Ncoii as shown in Fig 8(a). As seen in the plot, the Au+Au yield 
is higher than the one in d+Au in pj <2.5 GeV/c, which is close to what is expected in a literature [3]. When we 
divide Au+Au by A^co//-scaled d+Au data, we clearly see the existence of an additional effect in Au+Au collisions 
(Fig. 8(b)). 

Since both statistical and systematic uncertainties are still very large in d+Au measurement, it is hard to exhibit a 
concrete quantitative conclusion. However, in order to visually improve the result, we tried to parameterize the d+Au 
points by using a fit function. Fig 9(a) shows a fit to the d+Au points and Fig 9(b) shows the fit scaled by Ncoii and 
Au+Au points. The excess in Au+Au for pr <2.5GeV/c still looks to be definite compared to the yield in d+Au, 
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Figure 7: (a, left) Direct photon yield in d+Au collisions from this analysis together with Run3 real photon result and p+p yield scaled by Ncoii- (b, 
right) The same data with NLO pQCD calculation, instead of p+p yield. 
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Figure 8: (a, left) Direct photon yields in Au+Au and d+Au collisions scaled by the difference of Ncoii in both collision systems, (b, right) Ratio of 
Au+Au yield to d+Au yield scaled by Ncoii- 



T. Sakaguchi et al. /Nuclear Physics A 00 (2011) 1-8 



1 




(GeV/c) (GeV/c) 
Figure 9: (a, left) Fit to the d+Au invariant yield, and (b, right) comparison of the fit and Au+Au yield. 

implying that the excess is caused by a source created in hot dense medium. 

6. Toward better measurement of dileptons 

The statistical uncertainty and a part of systematic uncertainty in dilepton measurement are governed by the a 
background arising from random combinations of electrons from photons converted at beam pipe and Dalitz decay of 
jfi and 77. If we could tag and remove such electrons from foreground electron pair measurement, the errors would be 
reduced significantly. 

We developed a hadron blind detector (HBD) to realize this tagging. The detector is windowless Cerenkov type 
detector with CF4 gas in its radiator. We operate the detector in magnetic field free region in order to look for an 
opening angle of electron pairs. The principle of operation is following; if we see a cluster that has charge corre- 
sponding to the number of photo-electrons for single electrons, we keep the track associated with the cluster. If we 
see a cluster with double of the number of the photon-electrons, we reject the tract associated with the cluster, assum- 
ing that this is contributed by two electrons, which are likely by photon conversion or Dalitz decay. The detail of the 
detector description and its performance can be found in [7]. Fig. 10 shows the invariant mass distributions with and 
without using HBD information from Year-9 p+p run at V^=500GeV. The signal to background ratios in 0, cj and p 
meson region are significantly improved after applying a charge cut in HBD. We are planning to also apply the charge 
cut in HBD in Au+Au collisions in RHIC Year- 10 run data. It is essential in Au+Au collisions because the particle 
multiplicity is much higher than that in d+Au or p+p. 

7. Summary 

Electro-magnetic probes such as dileptons and photons are strong probes to investigate the thermodynamical state 
of the early stages of collisions since they leave the system unscathed. The PHENIX experiment has measured both 
photons and dileptons in p+p, d+Au and Au+Au collisions. An excess of dilepton yield over the expected hadronic 
contribution is seen in 0.2-0.8 GeV/c^ in Au+Au collisions, which is prominent in lower pj and most central. The 
invariant slopes are low for low mt - mo and high for high mt - mo, implying that the local slopes are roughly 
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Figure 10: Dilepton invariant mass spectra (a, left) without HBD matching and (b, right) with HBD matching and a charge cut. 



proportional to the energy of photons emitted. Direct photons are measured through their internal conversion to 
electron pairs. We saw a large enhancement in Au+Au collisions over p+p yield scaled by the number of binary 
collisions. It turned out from the latest results on d+Au collisions that this enhancement is not explainable by a 
nuclear effect. 
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